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ABSTRACT 

Non-radiating, advection-dominated, accretion flows are convectively unstable in the 
radial direction. We calculate the two-dimensional (r — 0) structure of such flows as- 
suming that (1) convection transports angular momentum inwards, opposite to normal 
viscosity and (2) viscous transport by other mechanisms (e.g., magnetic fields) is weak 
(a <C 1). Under such conditions convection dominates the dynamics of the accretion 
flow and leads to a steady state structure that is marginally stable to convection. We 
show that the marginally stable flow has a constant temperature and rotational velocity 
on spherical shells, a net flux of energy from small to large radii, zero net accretion rate, 
and a radial density profile of p oc r -1 / 2 , flatter than the p oc r -3 / 2 profile characteristic 
of spherical accretion flows. This solution accurately describes the full two-dimensional 
structure of recent axisymmetric numerical simulations of advection-dominated accre- 
tion flows. 



1. Introduction 



Analytical calculations of the structure of non-radiating, advection-dominated, accretion flows 
(ADAFs) have shown that they are convectively unstable in the radial direction (e.g., Begelman 
&i Meier 1982; Narayan & Yi 1994)0 In the absence of radiation, the entropy of the gas increases 
as it accretes (due to viscous dissipation) and the sub-Keplerian rotation of the flow is insufficient 
to stabilize this unstable entropy gradient. This conclusion has been confirmed by several numer- 
ical simulations (Igumenshchev, Chen, & Abramowicz 1996; Igumenshchev & Abramowicz 1999, 
hereafter IA; Stone, Pringle, & Begelman 2000, hereafter SPB). 

Narayan & Yi (1994; 1995) argued that convection was unlikely to significantly modify the 
structure of ADAFs, essentially because the inflow time of the gas would be shorter than the 
characteristic convective turnover time (i.e., advection would overwhelm convection). Their analysis 



1 Chandra Fellow 

2 Note that this is very different from thin accretion disks which can be vertically convective. 
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assumed that convection behaved like "normal" viscosity, transporting angular momentum and 
energy from small to large radii in the flow. 

A number of treatments of convection in thin accretion disks, however, have argued that 
convection does not act like normal viscosity; instead it transports angular momentum inwards, 
opposite to what is needed to allow accretion (Ryu & Goodman 1992; Stone & Balbus 1996; Balbus, 
Hawley, & Stone 1996). Narayan, Igumenshchev, & Abramowicz, (2000; hereafter NIA) have shown 
that, if present in ADAFs, this inward angular momentum transport by convection can dramatically 
modify the structure of the accretion flow, essentially suppressing accretion. In this paper we reach 
the same conclusion as NIA via a different analysis; we also find the full two-dimensional structure 
of the flow (NIA considered height integrated flows). 

We assume that convection transports angular momentum inwards and that outward transport 
by other mechanisms (e.g., magnetic fields) is weak in the sense that the induced radial velocities 
are highly subsonic (the Shakura-Sunyaev dimensionless viscosity or its magnetic analogue is small; 
a <C 1). In this case the steady state structure of the accretion flow must, statistically speaking, 
approach marginal stability to convection. This is because, in the absence of marginal stability, 
the characteristic convective turnover time is of order the dynamical time, much shorter than the 
inflow time of the gas (because a <C 1). Convection therefore dominates the flow structure and 
forces it to marginal stability (just like in the solar convection zone). 

In the next section (§2) we calculate the two-dimensional (r — 0) structure of an axisymmetric, 
radially self-similar accretion flow that is marginally stable to convection. This provides us with the 
angular structure of the flow. In §3 we argue that the only marginally stable accretion flow which 
satisfies mass, angular momentum, and energy conservation is one for which there is no accretion 
(M = 0), but a net flux of energy from small to large radii^j This is in contrast to standard ADAF 
solutions in which there is no net energy flux, but a finite accretion rate (see, e.g., Narayan & Yi 
1994; Blandford & Begelman 1999). In §4 we compare our solution to recent numerical simulations 
and discuss its relevance to "real" ADAFs. 



2. Marginal Stability: the Angular Structure 

We consider stability to local, axisymmetric, adiabatic, perturbations. Entropy gradients in 
ADAFs tend to destabilize perturbations while rotation stabilizes them. The net stability criterion 
can be found by requiring that the projection of the buoyancy and centrifugal forces on the dis- 
placement vector of a mode is negative (so there is a net restoring force for any perturbation); this 



3 In the self-similar limit, a finite energy flux is generated by a "black hole" of infinitesimal radius which accretes 
an infinitesimal amount of matter. 
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yields the "Hoiland" criterion (e.g., Tassoul 1978; Begelman & Meier 1982) 

(Va • dr)(g • dr) - ^ (V[i>*J2] • dr) d# < 0. (1) 

In equation (JlJ) , R = rsin# is the cylindrical radius, dr = drf + rdOO is the displacement vector, 
s = ln(p) — 7ln(p) is (7 — 1) times the entropy, g = —fv\/r + Rv^/R is the effective gravity, and 
vk = r -1 / 2 and are the Keplerian and rotational velocities, respectively (taking GM = 1, where 
M is the mass of the central object). 

Away from the boundaries, there is only one characteristic length scale in an ADAF, the 
distance from the central object (r). We therefore consider radially self-similar flows in which the 
dynamical variables are power laws in radius times functions of 9, i.e., p = r~ n p(6), c 2 = r^ 1 c 2 (6), 
p = r^ 1 ^ n p(9), and = r~ l v^{6); c 2 is the isothermal sound speed (i.e., temperature) andp = pc 2 
is the pressure. The radial scalings of the temperature and rotational velocity are set by the 1/r 
gravitational potential. For reasons that will become clear in §3 we allow the density to be an 
arbitrary power law in radius, p oc r~ n , rather than requiring the usual spherical flow value of 
n = 3/2. 

In the Appendix we show that, for a given n and 7, the requirement that the flow be marginally 
stable to convection for all 6 determines its structure.^ The temperature and rotational velocity 
are constant on spherical shells, i.e., independent of 9, with values 



C 2 = ^- and 4 = 2 V l[^\. (2) 



The density is then a power law in sin#, with 

p{6) oc (sin#) 2( ^- n) . (3) 

Note that in the above solution the surfaces of constant entropy and angular momentum coin- 
cide, and are given by rsin 2 # = constant. Such "gyrentropic" solutions are expected for marginal 
stability to convection, and have been discussed previously in the literature (e.g., Paczynski & 
Abramowicz 1982; Blandford, Jaroszynski, & Kumar 1985). 



3. Implications of Marginal Stability: the Radial Structure 



4 Marginal stability means that the quadratic form given by eq uatio n ([l]) is negative along all directions but one, 
along which it is equal to zero; see the discussion below equation (|A5[) in the Appendix. 
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The marginally stable solution derived in the previous section describes the angular structure 
of the flow for a given 7 and n. Here we argue that the requirement of marginal stability, together 
with mass, angular momentum, and energy conservation, also uniquely fixes the radial density 
power law, n. We first give a physical argument and then elaborate on it with some non-rigorous 
mathematics. 

Since radiation is assumed negligible, we write the conservation laws as 

dF f w 

— = 0, where F = dO sin 9F. (4) 
dr J 

F represents the time- aver aged radial flux of energy (Fe), angular momentum [F£), and mass (Fm) 
at a given r and 9, while F is the corresponding angle-integrated flux (i.e., the net flux through a 
spherical shell of radius r). 

In standard ADAF models, n = 3/2; as discussed below, this corresponds to a finite flux of 
mass onto the central object {Fm < 0), but zero net angular momentum and energy flux, i.e., 
Fe = Fl = (Narayan & Yi 1994; Blandford & Begelman 1999). That such a solution is possible 
can be understood as follows. Viscosity gives rise to a net torque on the gas which drives it 
inwards (Fm < 0). The outward angular momentum flux due to this torque is exactly balanced 
by the inward flux of angular momentum due to the bulk motion of the gas; this is how Fl = 0. 
Associated with the outward flux of angular momentum due to the viscous stress is an outward 
flux of energy due to the work done by the stress on the gas. The net energy flux is zero, i.e., 
Fe = 0, because the outward energy flux due to viscosity is exactly balanced by the inward flux of 
energy due to the bulk motion of the gas. The latter is given by FuBe, where Be is the Bernoulli 
constant, 

D 2 1 7 2 2 /r-\ 

2 7 — 1 

where v 2 = v 2 + v 2 + Vg . For standard ADAF solutions, Be must be positive and relatively large 
(i.e., a reasonable fraction of v\) in order for the inward mechanical flux of energy to balance the 
outward viscous flux (Narayan & Yi 1994; Blandford & Begelman 1999). 

From equation (||) one can readily show that the Bernoulli constant satisfies Be ~ for all 
9, 7 and n in the marginally stable state.^] It is therefore not possible to have a normal ADAF, i.e., 
an n = 3/2 solution. There is no inward mechanical flux of energy to balance the outward flux due 
to viscosity. 

Physically, the above argument shows that a marginally stable flow must have a net flux of 
energy to large radii. This uniquely fixes n for a given viscosity prescription since there is only 



5 An alternative proof that Be « is as follows: the marginally stable flow is gyrentropic, with surfaces of constant 
angular momentum and entropy coinciding. Under these conditions, one can show that the surfaces of constant 
Be coincide with the surfaces of constant angular momentum and entropy. Since v$ and c are independent of 9, 
Be is constant on spherical shells. Angular momentum and entropy are, however, constant on surfaces defined by 
r sin 2 6 = constant. Thus Be must be zero at marginal stability. 
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one n for which the outward flux of energy is independent of radius, as is required in steady state. 
Modeling viscosity via an r—cj) component of the stress tensor, the viscous energy flux is oc r 2 vpv^/r, 
where v is the kinematic viscosity. For a self-similar flow, the viscous flux is oc up oc ur~ n and 
so if v oc r@, a finite energy flux requires n = (3. For a Shakura-Sunyaev viscosity prescription, 
v = acr oc r 1 / 2 , so that n = 1/2. In this marginally stable solution with a finite energy flux, 
Fm = 0, i.e., there is no accretion; matter simply circulates in convective eddies. This is possible 
because the inward transport of angular momentum by convection balances the outward transport 
by viscosity. 

To elaborate on this discussion, we construct specific expressions for Fl and Fe- This elabo- 
ration is intended solely to motivate the various signs utilized in the above argument. We neither 
can nor do construct a fully satisfactory model for the energy and angular momentum fluxes in 
a turbulent/convective accretion flow (see NIA for an analysis using mixing length theory). The 
reader satisfied with the above argument may wish to proceed directly to §4. 

We consider a model in which the energy and angular momentum fluxes have contributions 
from two sources: (1) hydrodynamic effects such as bulk motion of the gas and convection and 
(2) non-hydrodynamic viscous processes such as those due to magnetic fields. We model the latter 
with an ^-component of the stress tensor, which we denote — T v , where T v > because "normal" 
viscosity transports angular momentum outwards. In the usual Shakura-Sunyaev prescription, 
T v = 1.5vpVd)/r, where v = acr is the kinematic viscosity. 

With this model, the radial fluxes can be written as 

F E = {r 2 pv r Be)+r 2 {v^)T v (6) 

F L = (pv r r 3 sin 6 V(j) ) + r 3 sin 9T V , (7) 

F M = (r 2 pv r ), (8) 

where ( ) denotes a time average over many convective turnover times. The first terms on the right 
hand side of equations @ and (0) are the hydrodynamic contribution to the energy and angular 
momentum flux; these will be explained in more detail below. The second terms on the right hand 
side of equations @ and (|^) are those due to the viscous stress T v . Associated with the angular 
momentum flux due to this stress is an energy flux oc QT V due to the rate at which the stress does 
work on the gas (fi = v^/fr s'mO] is the rotation rate of the gas). 

To understand the meaning of the hydrodynamic contribution to the angular momentum flux 
in equation (|7|), we decompose into mean and fluctuating components via = (v^) + v^. With 
this decomposition, 

(pv r r 3 sin Ovj,) = r sin 6(v^)F M + r 3 sin 6>T C , (9) 

where T c = (pVrV^) is the stress tensor associated with convection (the Reynold's stress); T c is 
negative because convection is assumed to transport angular momentum inwards. Physically, the 
two terms in equation (^) correspond to the angular momentum flux due to bulk motion and 
convection, respectively. 
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An analogous decomposition of the hydrodynamic contribution to the energy flux in equation 
(|6|) is less useful because the energy flux contains a large number of quantities that are third order 
in the fluctuating variables; a detailed analysis of these terms is non-trivial and not particularly 
illuminating. Instead, we simply write 

{r 2 pv r Be) = (Be)F M + r 2 {v^)T c + F c , (10) 

which is formally just a definition of the convective energy flux F c . Physically, the first two terms 
on the right hand side in equation (|To| ) correspond to the energy flux due to bulk motion of the 
gas and the energy flux due to the rate at which the "convective stress," T c , does work on the gas. 
The remaining term in equation (|l0|), F c , is the energy flux due to convection in the absence of 
bulk motion or angular momentum transport. This is, in fact, the energy flux usually considered 
in convective media (e.g., mixing length theory). In non-rotating stars, e.g., convection does not 
lead to a significant transport of mass or angular momentum, but it does lead to a significant flux 
of energy; in our notation this is captured by the term F c in equation ([n]). Although a detailed 
expression for F c in terms of the turbulent quantities is complicated, it is also unnecessary for our 
purposes. All we need is its sign; F c > because convection transports energy from high entropy 
(small radii) to low entropy (large radii). 

Using the Shakura-Sunyaev prescription for T v , one finds that the viscous contribution to 
the energy flux is oc r 1//2_n . Similarly the viscous contribution to the angular momentum flux is 
oc r 2 ~ n . In order to conserve energy, angular momentum, and mass (dF/dr = 0), the corresponding 
angle-integrated fluxes must either be zero or independent of radius. Therefore, we can write our 
final conservation equations as 

Fe = (Be)F M + F c + r 2 (^)T = Ci<5 n , 1/2 , (11) 

F L = r{v$)F M + r 3 TsW9 = C 2 5 n ,2, (12) 



Fm = {r 2 pv r ) = C 3 5 n:3/2 , (13) 

where the Ci are constants and T = T v +T C is the total stress tensor. In writing equations (|il"l)-([i"3l) 
we have integrated the fluxes on spherical shells and have used the fact that c 2 , v 2 ^, and Be are 
independent of 9 in the marginally stable state and thus can be taken out of the angular integration. 

Standard ADAF solutions correspond to n = 3/2. As indicated by equations (|1T1)-(|13|), this 
is a solution for which there is a finite flux of mass onto the central object, but no net energy or 
angular momentum flux (Narayan & Yi 1994; Blandford &: Begelman 1999). As argued above, 
no such solution is possible when convection dominates the structure of the flow and drives it to 
marginal stability. 



Consider a putative n = 3/2 solution to equations (|iT|)-(13). For such a solution Fm < since 
there is an inward flux of mass. From equation (|ll]), and using Be ~ near marginal stability, 
we see that Fe ~ r 2 v ( f ) T + F c . Both terms in this expression for Fe are positive and so it is 
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not possible to have Fe = 0-f\ There is consequently no solution with n = 3/2. This argument 
shows that an ADAF that is marginally stable to convection has a net flux of energy to large radii. 
Mathematically, the corresponding self-similar solution has n = 1/2. In this solution, Fyi = 0, 
i.e., there is no accretion; this is possible because the inward transport of angular momentum by 
convection exactly balances the outward transport by viscosity (T sin 9 = 0). 



4. Discussion 

Our argument can be summarized as follows. If convection in ADAFs transports angular 
momentum inwards and if outward transport due to other mechanisms is weak (a C 1), the time 
averaged structure of an ADAF must approach marginal stability to convection. By analyzing the 
structure of a marginally stable flow we find that there must be a net outward flux of energy due to 
convection. For a self-similar flow there can either be a net flux of mass or a net flux of energy, but 
not both (§3). Thus our marginally stable flow with a finite energy flux has zero accretion rate in the 
self-similar regime; the inward angular momentum transport due to convection cancels the outward 
transport due to viscosity. Physically, matter in a given spherical shell circulates indefinitely in 
convective eddies instead of accreting (in reality there will be a small net accretion rate due to the 
violation of self-similarity near the surface of the central object). This marginally stable solution 
has a radial density profile of p oc r -1 / 2 , rather than the usual p oc r~ 3 / 2 for spherical flows. 

SPB and IA have carried out axisymmetric hydrodynamical simulations of ADAFs; they find 
that the flow is convectively unstable for small a. Our analytical calculation reproduces the time 
averaged properties of these simulations rather well. SPB's run K corresponds to a Shakura-Sunyaev 
like viscosity with a = 10~ 3 ; in this simulation, they find p oc r -1 / 2 , as we derived in §3. SPB also 
consider non self-similar viscosity prescriptions for which v = p or v = constant. As discussed in §3, 
for v oc a finite energy flux requires n = [3 (where p oc r~ n ). This follows from the requirement 
that the viscous energy flux oc r 2 vpv^/r oc up be independent of radius. Thus for both v = p 
((3 = —n) and v = constant ((3 = 0), a net energy flux requires n = 0, i.e., p independent of radius. 
This is precisely what is found in the simulations (see SPB's Fig. 7). 

The interpretation of the numerical results as due to a marginally stable flow is supported by 
the excellent quantitative agreement between our results and the simulations.^ For example, SPB 
find rotational velocities that are nearly constant on spherical shells (see their Figs. 7 & 12), as 



An inward flux of mass with Fl — requires Tsin# > 0. In principle, one can have T < with Tsin# > since 
we do not know the 6 distribution of the convective (T c < 0) and viscous (T v > 0) stresses (this would correspond to 
an outward angular momentum flux and an inward energy flux due to "viscosity + convection"). It is likely, however, 
that Tsin^ > because T > at all 9, i.e., because viscosity wins out over convection (\T V \ > |T C |) at all 6. In this 
case T > follows. 

7 SPB noted that the angular structure of the contours of entropy, density, angular momentum, and other dynamical 
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we derived in §2 and the Appendix.^ For 7 = 1.5, our marginally stable n = 1/2 solution has 
c 2 = 0.22vjj-, in good agreement with IA's low a simulations, which find c 2 ~ Q.2v\ (see their Fig. 
10). In addition, our predicted angular density profiles of p oc sin 2 6 (for 7 = 5/3 and n = 1/2) and 
p oc sin 3 9 (for 7 = 5/3 and n = 0) are in excellent agreement with the simulations (see SPB's Fig. 
12 and Fig. 7, respectively), as is our prediction that the time averaged Bernoulli constant should 
be small, i.e., <C v 2 K (see SPB's Figs. 4 & 9 and the discussion before their eq. [9]). 

It is important to note that our analysis does not rely on the suggestion of Blandford &; 
Begelman (1999) that ADAFs will drive strong outflows. In fact, our density profile of p oc r -1 / 2 
is derived assuming no outflows. 

If relevant to real flows, the p oc r -1 / 2 density profile of the convection-dominated accretion 
flow would have important observational implications. Most importantly, it would naturally explain 
the low luminosities of many nearby supermassive black holes. For a given density and temperature 
in the ISM of a host galaxy, the accretion rate predicted by the convection-dominated accretion 
flow is ~ R s /Ra times smaller than the Bondi accretion rate, where R s is the Schwardschild radius 
of the black hole and Ra ~ GM/c\ is the accretion radius of the black hole (ca is the sound speed 
of the gas at large radii). This follows by evaluating M = AirR^pv close to the black hole where the 
velocity is of order the sound speed (the gas goes through a sonic point on its way into the black 
hole). Using p oc i? -1 / 2 and c oc i?" 1 / 2 then yields M ~ Mb(R s /Ra), where the Bondi accretion 
rate is Mb ~ ^R\paca (pa being the density of the gas in the vicinity of the accretion radius). 
Such a low accretion rate accounts naturally for the low luminosity of supermassive black holes 
such as Sgr A* at the center of our Galaxy. 

In addition, the X-ray to radio flux ratio and the X-ray spectrum of accreting supermassive 
black holes are sensitive to the variation of density with radius, i.e., n (Di Matteo et al. 2000; 
Quataert & Narayan 1999). Smaller n, e.g., n = 1/2 instead of n = 3/2, results in a larger X- 
ray to radio flux ratio and a much harder X-ray spectrum (dominated by bremsstrahlung instead 
of Comptonization) . This appears to be consistent with several low-luminosity AGN in nearby 
elliptical galaxies (Di Matteo et al. 2000). 

There are, however, several concerns about the relevance of the convection-dominated solution 
that can only be addressed by future numerical simulations: (1) Is the statistical steady state 
we have calculated unstable to non-axisymmetric perturbations? This is important because our 
conclusion that there must be a net energy flux follows from the properties of the marginally stable 



variables found in the simulations correspond to that expected near marginal stability (see the discussion in Begelman 
& Meier 1982). Moreover, Igumenschev, Chen, & Abramowicz (1996) noted that convection in low a flows should 
lead to a marginally stable state, which was consistent with the fact that the surfaces of constant specific angular 
momentum and entropy coincided in their simulations. 

8 Near the pole (8 ^ 10°), decreases in the simulations. This is because near the pole = constant corresponds 
to a divergent rotation rate. In reality, and in the simulations, this is smoothed out by a 8 — (f> component of the 
viscous stress tensor not considered in this work. 
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state. (2) Does radial convection transport angular momentum inwards in 3D flows? If not, one 
is likely to get a normal ADAF, perhaps modified by strong outflows (Blandford & Begelman 
1999; see also IA's high a simulations). Stone & Balbus (1996) have argued that outward angular 
momentum transport by convection requires coherent azimuthal pressure gradients. These are 
manifestly absent in the axisymmetric simulations of SPB and IA. Thus 2D convection is artificial 
in an important way. (3) How large is a due to magnetic fields? If a ^ 0.1, as is quite plausible, 
then convection may be a minor perturbation since the gas flows into the black hole on a timescale 
comparable to the convective turnover time. 

We are grateful to Ramesh Narayan for pointing out that inward transport of angular mo- 
mentum by convection was important for understanding the numerical simulations of Stone et al., 
and for several additional discussions. We thank John Bahcall for comments on the paper and 
the referee, Mitch Begelman, for several useful comments. AG was supported by the W. M. Keck 
Foundation and NSF PHY-9513835. EQ acknowledges support provided by NASA through Chan- 
dra Fellowship grant number PF9-10008 awarded by the Chandra X-ray Center, which is operated 
by the Smithsonian Astrophysical Observatory for NASA under contract NAS8-39073. 
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A. Calculation of the Marginally Stable Solution 

In this Appendix we calculate the angular structure of a marginally stable flow by analyzing 
equation (|l|) assuming radial self-similarity. Equation (|l|) can be written as the inequality 

Ax 2 + By 2 + Cxy < 0, (Al) 

where x = dr and y = rdd. 

The coefficients of the quadratic form are given by 

A = g r d r s - ^ sm9d r (R V4> ), (A2) 

B = ^d e s - ^ coe0fy(ifov), (A3) 
r rH £ 

C = g e d r s + ^dgs - ^ (cosOdriRvj,) + ^-OeiRv^ . (A4) 

The marginally stable flow corresponds to 

A < 0, B < 0, and C 2 = AAB. (A5) 

The first two conditions follow because the flow must be stable to x = or y = perturbations 
(purely tangential and purely radial, respectively). The condition C 2 = 4AB follows because the 



and 



left hand side of equation (Al) defines a quadratic function with negative curvature. If it has 
two real roots then it is positive somewhere and hence there are perturbations which are unstable. 
Marginal stability occurs when the two roots collapse to one (C 2 = 4^45) and stability happens 
when the two roots of the quadratic function are imaginary (C 2 < <iAB). 

In addition to marginal stability, the flow must satisfy radial and 9 momentum balance 

i = - p -% p + i 

r r 
and 

cot 9pvl = d e p, (A7) 
where p = pc 2 is the pressure and c is the isothermal sound speed. 



Vk - p- 1 d rP + V -* (A6) 



With radial self-similarity, equations (A6) and ( |A7| ) become 

1 = (1 + n)c 2 + v\ (A8) 

and 

cot 9 pvl = c 2 p' + p(c 2 )', (A9) 
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where ' = dg and all dynamical variables (e.g., p, c 2 ) are now dimensionless functions of 9. 
For a self-similar flow the coefficients of the quadratic form become 

A = 1 — 71(7— 1) + v l («(7 " 1) " (7 + 1)) , (A10) 

B = ((1 " 7) cot fc;J(l + n) - 2 7 v^) - |^ cos 2 + ^ cos sin 0) , (All) 

and 

C = 2uJcat0((7-l)n-(7 + l)). (A12) 

The marginal stability requirement is that C 2 = AAB. Defining a new independent variable, 
£ = — In sin 0, this can be written as an equation describing the marginally stable temperature 
profile 

dc 2 (c 2 -c 2 )(c 2 -c 2 ) 

dc - 4^ ' (A13) 

where 

c\ = — ^— , 4 = —, c l = c A±A. (A14) 

7-1 

Since the region extending from the equator to the pole corresponds to £ going from to 00, the 
only non-singular solution which extends from the equator to the pole has c 2 = c 2 . Any other 
solution becomes singular somewhere between the equator and the pole or has a free surface at a 
finite (see below). 

The above solution has c 2 and constant on spherical shells, i.e., independent of 6, with 
values 

1 i — fi 



From equation (A7) the density is then a power law in sin#, with 

poc (sin#) 2( ^~ n) . (A16) 



Equation ( A15| ) gives the marginally stable state corresponding to C 2 = 4AB. One can show 



that this state also satisfies A < and B < and so all of the conditions in equation (A5) are 
satisfied. 

The referee has pointed out, correctly, that there exists another class of marginally stable 
solutions not considered in this paper. These have a free surface at a finite 6, where the density 
and pressure vanish. Outside of the free surface lies an empty funnel, where a hydrodynamic or 
magnetohydrodynamic outflow may be driven. The properties of these solutions are discussed in 
Blandford & Begelman (in preparation). 
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One final comment on convective stability is in order. In evaluating the stability of height 
integrated ADAFs to convection, Narayan & Yi (1994) and NIA used iV 2 + k 2 > as their stability 
requirement, where N is the Brunt- Vaisala frequency and k is the epicyclic frequency. In our 
notation, this is equivalent to taking A < since, at the equator, A = — 7(iV 2 + k 2 ). In the full 
2D flow, the criterion A < is, however, a necessary, but not a sufficient, criterion for convective 
stability (see the discussion below eq. flAlJ). Physically, A < is the stability requirement for 
radial perturbations, i.e., for y oc dO = (see eq. flAll ). Near marginal stability, however, the 
"almost" unstable modes have y = —^JAjBx and are thus not radial perturbations. True marginal 
stability therefore corresponds to C 2 = AAB and not A = 0. It turns out that our marginally flow 
has a rotational velocity a factor of v2 larger than the flow defined by A = 0.0 NIA found that 
their marginally stable flow (defined by A = 0) underpredicted the rotational velocities found in 
IA's simulations by a factor of y/2 (see the discussion in their §7); this is readily explained by 
the above analysis, 



This additional rotation is needed to stabilize all perturbations. 

10 To be concrete, for 7 = 5/3, our n = 1/2 solution has c 2 = 0.28«x, while NIA's solution gives c 2 = 0A8v^; the 
simulations find c 2 ~ 0.25t>|-, in excellent agreement with our result (see §7 and Fig. 4 of NIA). 



